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Nonlinear effects in superconductors are important for both basic science 
and technology. A strong nonlinearity may be exploited to demonstrate some 
important quantum phenomena in the microwave region, such as quantum 
squeezing 1,2,3 and experimental observation of the so called dynamical Casimir 
effect 4 ; Whereas technologically, these effects may allow some intriguing applica- 
tions such as bifurcation amplifiers for quantum measurements 5,6 and resonant 
readout of qubits 7 . In this work we study the response of a superconducting 
(SC) microwave stripline resonator to a monochromatic injected signal. We find 
that there is a certain range of driving parameters, in which a novel nonlin- 
ear phenomenon immerges, and self-sustained modulation (SM) of the reflected 
power off the resonator is generated by the resonator. That is, the resonator 
undergoes limited cycle oscillations, ranging between several to tens of mega- 
hertz. A theoretical model which attributes the SM to a thermal instability 
yields a good agreement with the experimental results. A similar phenomenon 
was briefly reported in the 60 's 8, &i&ii in dielectric resonators which were par- 
tially coated by a SC film, but it was not thoroughly investigated and therefore 
its significance was somewhat overlooked. This phenomenon is of a significant 
importance as it introduces an extreme nonlinear mechanism, which is by far 
stronger than any other nonlinearity observed before in SC resonators^ 2 -. It 
results in high intermodulation gain, substantial noise squeezing, period dou- 
bling of various orders 12 , and strong coupling between resonance modes (see 
supplementary figure 1). 

Our device (Fig. [I] (b) — (d)) integrates a narrow microbridge into a SC stripline ring 
resonator. The impedance of the microbridge strongly affects the normal modes of the res- 
onator and thus the resonance frequencies can be tuned by either internal (Joule self heating) 
or external (infrared illumination^ 3 -) perturbations 1 -. Further design considerations, fabrica- 
tion details as well as normal modes calculation can be found elsewhere^ 3 -. The experiments 
are performed using the setup described in Fig. HJa). We stimulate the resonator with a 
monochromatic tone at an angular frequency u p , and measure the reflected power off the 
resonator by a spectrum analyzer in the frequency domain and an oscilloscope in the time 
domain. Measurements are carried out while the device is fully immersed in liquid Helium. 

Fig. El panels (a) and (b) show typical experimental results of the SM phenomenon in the 
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FIG. 1: (a) SM measurement setup, (b) Schematic layout of the device. The resonator is designed 
as a strip line ring, having a characteristic impedance of 50 £1. It is composed of 200nm thick NbN 
deposited on a Sapphire wafer. A weakly coupled feedline is employed for delivering the input and 
output signals. An optical microscope image of the ring resonator section at which the microbridge 
is integrated is seen in panel (c), whereas panel (d) shows the microbridge, whose dimensions are 
1 x 10/mr 2 . 



frequency domain. The resonator is stimulated with a monochromatic tone at the resonance 
frequency of the third mode fs, and the dependence of the SM on the injected pump power is 
investigated. At relatively low and relatively high input power ranges (P pU mp ^ —33.25 dBm 
U -Ppump > —23.7 dBm) the response of the resonator is linear, namely, the reflected power off 
the resonator contains a single spectral component at the frequency of the stimulating pump 
tone uj p . In between these power ranges, regular SM of the reflected power off the resonator 
occurs (see panel (b), subplots (ii, iv)). It is realized by rather strong and sharp sidebands, 
which extend over several hundred megahertz at both sides of the resonance frequency. 
The SM frequency, which is defined as the frequency difference between the pump and the 
primary sideband, increases as the pump power increases. The regular SM starts and ends 
at two power thresholds, referred to as the lower and the upper power thresholds. The lower 
power threshold (panel (b), i) spreads over a very narrow power range of approximately 
10 nW, during which the resonator experiences a strong amplification of the noise floor 
(noise rise) over a rather large frequency band, especially around the resonance frequency 
itself. This noise rise can be explained in terms of nonlinear dynamics theory, as it predicts 
the occurrence of strong noise amplification near a threshold of instability^^. The upper 
power threshold (panel (b), vi) spreads over a slightly larger power range than the lower one 
and has similar, but less extreme characteristics. 

The SM nonlinearity is also strongly dependent on the pump frequency. This dependence 
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FIG. 2: Typical experimental results of the SM phenomenon in the frequency domain. Panel 
(a) plots a color map of the reflected power P re fl as a function of the pump power P pump and the 
measured frequency /sa ; centralized on the pump frequency, which coincides with the the resonance 
frequency = 5.74GHz (/g A = /sa — /a)- Panel (b), subplots (i) , (ii) , (iv) , (vi) , plot the same 
measurement at four different pump powers, corresponding to (i) the lower power threshold, (ii) 
and (iv) powers that are in the range of regular SM, and (vi) the upper power threshold. The red 
curve in subplot (ii), which has a Lorenzian shape, shows the spectral density of the SM sidebands 
as predicted theoretically by the model (Eq. 26 of Ref. 15 ). The solid red curve in subplot (i) 
was obtained by numerically integrating the equations of motion of the model with a nonvanishing 
noise at the first power threshold and evaluating the spectral density^. Subplots (iii) and (v) were 
obtained by numerically integrating the equations of motion of the model for the noiseless case, 
at input powers corresponding to subplots (ii) and (iv) respectively, and calculating the spectral 
density. 

is shown in the supplementary video 1, where each frame is a graph similar to Fig. [2l 
panel(a), and corresponds to a different pump frequency, starting at a red shifted pump 
frequency (Af = f pump — fa < 0) and ending at a blue shifted pump frequency (Af > 0). 
Fig. [3] plots the SM frequency as a function of the pump frequency and power. The inset 
plots the stability diagram of the resonator (to be discussed below) as a function of the 
same parameters. The SM occurs only within a well-defined frequency range around the 
resonance frequency. A small change in the pump frequency can abruptly ignite or quench 
the SM. Once started though, the modulation frequency has a weak dependence on the 



FIG. 3: SM frequency / sm as a function of the pump power and the pump frequency, centralized on 
the resonance frequency /p Ump = /pump — /3 • The inset shows a schematic diagram of the stability 
zones of the resonator as a function of the pump power and frequency. The green, pale blue, yellow, 
and grey colors represent the SC monostable (MS(S)), NC monostable (MS(N)), bistable (biS), 
and astable (aS) zones respectively. 



pump frequency. The maximum SM frequency measured with this device is approximately 
41 MHz, whereas higher maximum of approximately 57 MHz has been measured with other 
devices^. 

We propose a theoretical model according to which the SM originates from a thermal 
instability in the microbridge section of the SC stripline resonator. Current-carrying su- 
perconductors are known to have two or more metastable phases sustained by Joule self- 
heating^. One is the SC phase and the other is an electrothermal local phase, known as 
hotspot, which is basically an island of normal conducting (NC) domain, with a temperature 
above the critical one, surrounded by a SC domain. This phenomenon can be explained by 
the heat balance equation holding at more than one temperature. A perturbation can trig- 
ger or suppress the formation of a hotspot and thus the microbridge can oscillate between 
instable phases. Such oscillations were often observed in experiments, for the case of a SC 
microbridge driven by an external dc voltage (see review-^ and references therein). 

In the current the microbridge is integrated into a stripline resonator, the system 

is driven into instability via externally injected microwave pump tone. We herein briefly 
present the corresponding theoretical model, whereas the full derivation of the equations, as 



well as experimental justifications to several assumptions are included elsewhere^ 5 -. Consider 
a resonator driven by a weakly coupled feedline carrying an incident coherent tone b m e~' lWpt , 

I ■ 1 2 

where b m is constant complex amplitude (\b m \ oc P pU mp) and u p is the driving angular 
frequency. The mode amplitude inside the resonator can be written as Be~ luJpt , where B (t) 
is a complex amplitude which is assumed to vary slowly on a time scale of l/u p . In this 
approximation, the equation of motion of B reads 2 . 



where ujq (T) is the temperature dependant angular resonance frequency, T is the tempera- 
ture of the hotspot, 7 = 71 + 72, where ji is the coupling constant between the resonator and 
the feedline, and 72 (T) is the temperature dependant damping rate of the mode. The term 
c m represents an input Gaussian noise with a zero-mean and a random phase. We consider a 
case where the nonlinearity originates by a local hotspot in the microbridge. If the hotspot 
is assumed to be sufficiently small, its temperature T can be considered as homogeneous. 
The temperature of other parts of the resonator is assumed to be equal to that of the coolant 
T . The power Q heating up the hotspot is given by Q = nQ t where Q t = tko 2'y2 \B\ is 
the total power dissipated in the resonator, and ^ k ^ 1 represents the portion of the 
dissipated power that is being absorbed by the microbridge. The heat balance equation 
reads 



where C is the thermal heat capacity, W = H (T — T ) is the heat transfer power to the 
coolant, and H is the heat transfer coefficient. 

Nonlinearity, according to our simple theoretical model, results by coupling the equation 
of motion of the mode amplitude in the resonator, Eq. (CD), to the thermal balance equation 
in the microbridge, Eq. (j2J) . The coupling mechanism is based on the dependence of both the 
resonance frequency ujq and the damping rate 72 of the driven mode on the resistance of the 
microbridge 14 , which in turn depends on the temperature of the microbridge. We assume the 
simplest case, where this dependence is a step function that occurs at the critical temperature 
T c . We further assume that all other parameters are temperature independent. In our 
previous publication 13 we have thoroughly investigated the dependence of the resonance 
modes on the resistance of the microbridge. We have shown that a resonance frequency can 
be smoothly shifted by gradually controlling the bridge resistance. In addition the damping 
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rate increases up to a certain maximum, as the resistance increases, but then decreases back 
to approximately its original value, when the resistance is further increased. 

Due to the step-like dependence of the resistance, and the dependence of the heat gen- 
eration, on the temperature of the microbridge^, the system may have in general up to 
two locally stable steady states, corresponding to the SC and NC phases. These states are 
associated with stationary solutions of the uncoupled Eq. (1) and (2), where the noise is 
disregarded. The stability of each of these phases depends on both the power and frequency 
parameters of the injected pump tone. In general there exist four different stability zones 
(see inset in Fig[3]) 15 . Two are monostable zones, where either the SC phase or the NC phase 
is locally stable. Another is a bistable zone, where both phases are locally stable^ 1 ^. The 
third is an astable zone, where none of the phases are locally stable, and consequently the 
resonator is expected to oscillate between the two phases. As the two phases significantly 
differ in their reflection coefficients, the oscillations translate into a modulation of the re- 
flected pump tone. The mechanism introduced in this theoretical model is somewhat similar 
to one of the mechanisms which cause self oscillations in an optical parametric oscillator—. 

The dependence of the SM on the input power, as observed in the time domain (Fig. 
SJ), gives a complimentary understanding of the SM phenomenon. Below the lower power 
threshold the resonator is in the SC low reflective phase, where only a small portion of 
the injected power is reflected off the resonator. Once the pump power approaches the 
power threshold, sporadic spikes in the reflected pump power occur, as seen in subplot (i). 
Near the threshold the device is in sub-critical conditions and these spikes are caused by 
a stochastic noise which triggers transitions of the microbridge from the SC phase to the 
NC high reflective phase. During a transition the stored energy in the resonator is quickly 
discharged and dissipated. Consequently, after a quick and short worm-up the microbridge 
cools down until eventually it switches back to the SC phase. As a result, the damping rate 
decreases back to its original value and thus the stored energy in the resonator is slowly built 
up again. Once a spike is triggered, the noise has a negligible effect on the dynamics of the 
energy discharge and buildup and therefore the line-shapes of the various spikes are similar. 
This behavior can also be clearly observed in the simulation results, shown in subplot (v). 

Regular SM of the reflected power (Fig. H] (ii,vi)) occurs when the pump power is set 
above the lower power threshold. In this case the pump tone drives the resonator to the 
astable state and the noise has a negligible influence. The dynamics of the oscillations is 




FIG. 4: Panels (i) — (iv) show typical experimental results of the SM phenomenon in the time 
domain. The reflected power is normalized by the value of the maximum peak to peak amplitude. 
AC coupling is employed in the measurement and thus a zero value represents the average measured 
power. Each panel refers to a different pump power range, (£) lower power threshold, (ii) and 
(Hi) powers corresponding to regular SM. (iv) upper power threshold. Panels (v) — (viii) show 
numerical integration of Eqs. (1) and (2) calculated for the above-mentioned four cases. The 
graphs are centralized around their mean value and normalized by the value of the maximum peak 
to peak amplitude. In the simulation of the lower and upper thresholds (panels (v) and (viii)), a 
Gaussian thermal noise, corresponding to 4.2 K and 15 K temperatures was assumed respectively 
, whereas in the simulation of the regular SM, (panels (vi) and (vii)), the noise was disregarded. 
Supplementary video 2 shows the SM in the time domain while gradually increasing the pump 
power from frame to frame. Each frame has a subplot similar to panels (i)— (iv), and a subplot 
similar to[2j6) panels (ii), (iv). 



similar to the one just described where the energy buildup time is strongly dependant on 
the pump power. Thus, the oscillation frequency is faster for higher pump power values 
(Fig. H] (Hi, vii)). The upper power threshold (Fig. H] (vi)) resembles the lower one, but 
the SC and NC phases exchange roles. The resonator is in the NC high reflective phase, 
and noise-induced spikes temporarily drive it to the SC low reflective phase. The internal 
thermal noise in the upper threshold is stronger than in the lower one and consequently this 
power threshold range is wider. 



In conclusion, we report on a novel nonlinear phenomenon where SM is generated in a SC 
microwave stripline resonator. This behavior is robust and occurs in all of our devices and at 
various resonance frequencies. A theoretical model according to which the SM originates by 
a thermal instability is introduced. In spite of its simplicity, the model exhibits a good quan- 
titative agreement with the experimental results. In our measurements we find strong and 
correlated nonlinear phenomena accompanying the SM^. One of these phenomena is self- 
excitation (SE) of coupled resonance modes, which is shortly described in a supplementary 
figure. This phenomenon may be exploited for state readout in quantum data processing 
systems 3 -. 
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